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Abstract 

We calculate the transition density for the overtone of the isoscalar giant monopole resonance 
(ISGMR) from the response to an appropriate external field ^ /g(r) obtained using the semiclas- 
sical fluid dynamic approximation and the Hartree-Fock (HF) based random phase approximation 
(RPA). We determine the mixing parameter ^ by maximizing the ratio of the energy- weighted sum 
for the overtone mode to the total energy-weighted sum rule and derive a simple expression for the 
macroscopic transition density associated with the overtone mode. This macroscopic transition 
density agrees well with that obtained from the HF-RPA calculations. We also point out that the 
ISGMR and its overtone can be clearly identified by considering the response to the electromagnetic 
external field ~ joiqr). 

PACS numbers: PACS numbers : 21.60Jz, 24.30.Cz, 26.60.Ev, 24.10.Nz 
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I. INTRODUCTION 



The properties of a giant resonance in nuclei are commonly determined from the distorted 
wave Born approximation (DWBA) analysis of its excitation cross-section by inelastic scat- 
tering of a certain projectile. The transition potential required in actual implementation of 
DWBA calculation is usually obtained by convoluting the projectile-nucleus interaction with 
the transition density associated with the giant resonance. The relevant transition density 
can be obtained from a microscopic theory of the giant resonance, such as the Hartree-Fock 
(HF) based random phase approximation (RPA). However, the use of a macroscopic tran- 
sition density p^'^^lr) greatly simplifies the application of the giant multipole resonance 
theory to the analysis of the experimental data. The simple form of the transition density 

Ptr(r)=«o(3 + r|-^Peq(r)roo(f), (1) 

obtained by the scaling approximation, is a well-known example of the macroscopic transi- 
tion density p™'*'^''(r) commonly used in the case of the isoscalar giant monopole resonance 
(ISGMR) l|. The transition density of Eq. nicely agrees with the ISGMR transition 
density obtained in microscopic HF-RPA calculations. It has a one-node structure, satisfy- 
ing the condition of particle number conservation J pl'^{r) dr = 0. Unfortunately the scaling 
consideration can not be extended to the overtone of the ISGMR, where p™'^'^''(r) has a two- 
node structure. To derive the macroscopic transition density p™^'^''(r) in this more general 
case, one can use the well-known method of determining p™'^'^''(r) from the local sum 

rule which is exhausted by one collective state with the appropriate choice of the transition 
operator /^(r). However, in the quantum random phase approximation, the highly excited 
collective modes are strongly fragmented over a wide range of energy and a special averaging 
procedure must be employed to determine the macroscopic transition density corresponding 
to an average collective excitation. In this respect, the semi classical Fermi- liquid approach 
(FLA) P] is more appropriate. Both the main ISGMR and its overtone are well-defined 
within the FLA as single resonance states. This fact enables us to derive the transition op- 
erator /^(r) simply by maximizing the fraction of the energy- weighted sum rule (FEWSR) 
exhausted by the single overtone. 

In the present work we suggest a procedure to derive the macroscopic transition density 
for the ISGMR overtone using both the HF based RPA and the Fermi-liquid approaches. 



We remark that some preliminary results of this investigation were presented in Ref. [5| (see 
also Ref. ^). 



II. QUANTUM DERIVATIONS 



The transition density ptr,n(r) for a certain eigenstate \n) of a nucleus with A nucleons 
is given by 

Ptr,n(r) = (0|p(r)|n), (2) 

where p{v) = Ylf=i '^(^ "~ is the particle density operator and |0) represents the ground 
state of the nucleus. The transition density reflects the internal structure of the nucleus and 
does not depend on the external field. However, a problem arises if one intends to derive the 
transition density ptr(i') for a group of the thin-structure resonances in the giant multipole 
resonance (GMR) region. An appropriate averaging procedure is necessary in this case and 
Ptr(r) can be evaluated if the nucleus is placed in an external field 

A 

Kxt- F(M) = ^/(r,), (3) 

i=l 

where the transition operator F({rj}) is so chosen that it provides a preferable excitation 
of the above mentioned thin-structure resonances. 
Let us introduce the local strength function 

5(r,E) = J] (0|p(r)|n) {n\F\Q) 5{E - E^) (4) 

n 

and the energy smeared local strength function S'(r, E) defined near the GMR energy Eji 
by 

_ 1 j-ER+AE/2 

Siv,E^) = —; dES{v,E). (5) 

JEr-AE/2 

The corresponding strength functions are given by 



S{E) = Jdr /(r) 5(r, E) = J2 J dr />) Ptr,„(r 



6iE-E^), (6) 



and 

S{E)= I dr /(r) S{r,E). (7) 



Let us assume, for the moment, that the operator F({rj}) excites only a single state \D), 
within the energy interval ± AE/2. The corresponding transition density ptr,D(r) = 
(0|p(r)|D) is then given by the following exact expression 

AE 



Ptr,Dlrj 



--S{r,En). 



(8) 



'S{En) AE 

We will extend expression (jHJ to the case of a group of the thin-structure resonances in the 
GMR region which are excited by the operator F({rj}) and define the smeared transition 
density ptr,R(r) as 

Ptr,R(r) = , ^(r. En). (9) 

^JS{ER) AE 

Note that Eq. Q is associated with the strength in the region of Er^AE /2 and is consistent 
with the smeared strength function S{Er) for a single resonance state. That is (see also Eq. 



S{E, 



1 



AE 



dr /(r) ptr,R(r) 



We also point out that with the Lorentz's function 
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g'fiE,En) 



the energy smeared S{r, Er) is given by 



Sir, E 



R) 



n{E-ERy + ^^' 



dE SiT,E) g^{E,ER) 



(10) 



(11) 



and the smeared transition density ptr,R(i') is obtained from 
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Ptr,R(r) = — 

\JS{Er) Tcy 

The consistency condition, Eq. (jiup . then reads 



--S{r,ER). 



(12) 



(13) 



S{Er) = — 
7r7 



dr /(r) Ptr,R(r) 



(14) 



In the quantum RPA, the local strength function 5'(r, E) is related to the RPA Green's 
function G{r',r,E) by 



Sir,E)= / />' 



-ImG(r',r,E) 
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dr' . 



(15) 



For the isoscalar monopole and dipole excitations, the transition operator /(r) is taken in 
the form of 

/(r) = /^(r) = /5(r)yoo(f) for L = 0, (16) 

and 

/(r) = /,(r) = /,(r)rio(f) for L = 1, (17) 

with an appropriate choice of the radial functions f^{r) and fr^{r), see below. In the following, 
the quantum transition density for the main ISGMR and its overtone is evaluated using the 
Eq. with Eft ± Ai?/2 taken separately for the ISGMR and the overtone regions. 

III. MACROSCOPIC TRANSITION DENSITY 

Let us consider the local energy- weighted sum Mi(r) given by (see Eq. (jH)) 



Mi(r) 



poo 

/ dE ES{r,E) = V E„, (0|/5(r)|n) (n|F|0), (Eq = 0). (18) 

-'0 



The continuity equation provides the following sum rule 



M,M = ~ (ftc,V/(r)) . (19) 

Let us assume that only one state \D) exhausts the sum rule Eq. ()19|). Then for the 
corresponding ("macroscopic") transition density, p'^^'^^{r), we have from Eqs. (fTHj) and ()19p 
the following expression 

pr^(r) = {0\p{r)\D) = ao V (peqV/(r)) , (20) 

where the normalization coefficient ao can be found using the energy-weighted sum rule 
(EWSR) 

m, = jdv /(r) Mi(r) = ^ j dv peq (v/(r))' . (21) 



Taking into account that 



{D\F\0) = dv mp^'ir) 



we obtain 

|ao| = ih^/2m){EDmi)-^^^. (22) 



Thus, the macroscopic transition density p™'^"(r) of Eq. ()20|) coincides with the quantum 
transition density for a certain state \D) if the single state \D) exhausts all the EWSR, Eq. 
(|2ip. associated with the transition operator /(r). In the case of the transition operator 
/(r) from Eqs. and (jl7|) . the EWSR is given by 







1) ( - 







(23) 



Assuming that the energy-weig hted transition strength (EWTS) Ed \{D\F\0)\^ for the 
resonance at Ed fully exhausts the EWSR associated with /(r) = f{r)YLM{^), we obtain 
the expression for the macroscopic transition density for the state at Ed from Eq. ()20j) 

L{L + 1] 



PtT"(r) 



2m 



2L + 1 



r dr 



Peq(r)yLM(f). (24) 



TTIiEd 

For L = we have from Eq. ()24|) the commonly used ISGMR result of Eq. (^. The main 
(lowest) ISGDR is the spurious state with the eigenenergy Ei- = 0. The next ISGDR is the 
overtone. Assuming that the EWTS for the 1~ overtone equals the EWSR associated with 
the operator 

/,(r) = (r3-r/r)rio(i-), (25) 
we obtain from Eq. (j2Up the macroscopic transition density as 



Pt: 



r .overtone \ 



V3 (^10 r + (3 - 77)^^ peq(r)Yio(f) 



for 



(26) 



For the isoscalar dipole mode, the translation invariance condition is used for the derivation 
of 77. This condition implies that the center of mass of the system can not be affected by 
internal excitation. We thus have, 



dr r F;o(f-)Pt 







for 



1. 



(27) 
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3, 



(28) 



From Eqs. (|^ and one obtains, see also Refs. 

where (r^) = dr r^peq(r)/ dr r'^peq{r) is the mean square radius. 

The free (mixing) parameter appearing in the transition operator /g(r) (similar to rj of 
Eq. (fT7|) for the L = 1 case) can be determined by an appropriate condition leading to a 
general method for the evaluation of the transition density for the overtone mode. In this 



work we present this method for the case of the monopole mode, L = 0. Let us introduce 
the transition operator /^(r) as 

/e(r) = (/-er')yoo(f) for L = 0. (29) 

The corresponding macroscopic transition density p'^^'^^{r) is obtained from Eqs. (j^UI) and 
(Uni) as 

Pcq(^)^oo(f) for L = (overtone). 

(30) 

The determination of the parameter ^ in Eq. requires an additional consideration since 
for the L = case we have no fundamental condition such as Eq. (jTTj) for the L = 1 mode. 
We note, however, that if we assume that the ISGMR has the transition density of Eq. ([T} 
and require that 

I d4(r)p-(r) = 0, (31) 
i.e., the ISGMR is not excited by the scattering operator of Eq. we have 

e = 2{r')/{r'). (32) 

Similar result is obtained by imposing the condition that the scattering operator r^Yoo{r) 
does not excite the overtone of the ISGMR, assuming the transition density of Eq. (jHOl)- 

Following the general requirement for the proper use of Eq. (pUj) in the derivation of 
the macroscopic transition density p™'^'^''(r), we can determine the parameter ^ from the 
condition that the transition operator /^(r) provides for the single overtone the maximum 
fraction of the energy- weighted sum rule mi of Eq. (|2T|l . 



macr /„\ 

Ptr (r) 



2aci 



lOr^ - 3^ + r (2r2 - ^) 



dr 



IV. SEMICLASSICAL FERMI-LIQUID APPROACH 

The transition density Ptr(r) and the strength function S{E) can also be evaluated within 
the semiclassical Fermi-liquid approach. For a given multipolarity L and overtone n, the 
FLA transition density is given by |^ 

Pfr'fn(r) = "Ln [^(i?o " r)ji(gi„r) + 



7 



Pol^Lo(f), (33) 



qLr. 

where po is the bulk density, Rq is the equihbrium nuclear radius and the parameter a is 
determined by the translation invariance condition (Eq. (j2ZI)) in the case of the isoscalar 
dipole compression mode and is given by 

a = ji{x)/xj[{x), X = qmRo- (34) 

The wave numbers qin are derived from the boundary conditions of the FLA model: the 
normal component of the tensor pressure 6P^f^, created by a sound wave, on the free surface 
of the nucleus should be equal to zero 

5Prr\n=R, = 0. (35) 

Note that for the case of compression sound modes, the contribution from the surface tension 
pressure is negligible and it was omitted in Eq. (IH^ . The boundary condition (Eq. 
leads to the following secular equation (see Ref. jl5l |) 

4 u 

[qr joiqr) - D^ji{qr)]r=Ro = 0, D^ = ^, (36) 

mpoCo 

where cq is the zero sound velocity and the coefficient p determines the contribution from the 
dynamical Fermi surface distortion associated with the collective motion in a Fermi liquid. 
In the case of a quadrupole distortion of the Fermi surface, one has Q| 

/. = Im (— ^)P,q. (37) 

1 — tUJT 

Here, Pgq ~ (2/5)ei?po is the equilibrium pressure of a Fermi gas, ep is the Fermi energy, and 
uj and T are the eigenfrequency and the relaxation time for sound excitations in the Fermi 
liquid, respectively. The relaxation time r is assumed to be frequency dependent because of 



the memory effect in the collision integral 



IS assumea to De irequency a 
]. Following Refs. \m,\m 



19| we take 



r = 47r'f3h/{hujy, (38) 

where f3 is the constant related to the differential cross section for the scattering of two 
nucleons in the nuclear interior. In the case of isoscalar sound mode, we will adopt (3 = 1.5 
MeV . The eigenfrequency u is obtained from the dispersion relation 

uj^ - clq^ + icu7g^ = 0, (39) 



where cq is given by 



cl = ^{K+12fi/po). (40) 

Here K is the nuclear incompressibihty coefficient and 7 is the friction coefficient 

4z/ ^ T \ 
7 = o ' iy = Re{ — Peq. (41) 



3pom \1 ~ '''^'^ 

The smeared FLA strength function S^^^{E) can be obtained in a way similar to S{E), 
obtained within the quantum approach of Eq. (fT^ . That is, 

2 



9,{E,ELn). (42) 



The smearing function g^{E, E^n) in Eq- is given by (72,^ <C -^Ln) 

Here, 7l„ is the damping parameter due to the viscosity of the Fermi liquid and Ein = 
hRe{ujLn), where the eigenfrequency Un is obtained as a solution to both the dispersion 
equation (j39|) and the secular equation (j36|l . We point out that the amplitude am in Eq. 
dnSl) for the FLA transition density Pt^Lni''^) is derived as the amplitude of the quantum 
oscillations 

am = ^/n/2BL{q) uj{q), (44) 

where q = qm is determined by Eq. (jHBjl and Bi^q) is the corresponding mass coefficient 
with respect to the density oscillations. The collective mass coefficient BL^q) can be found 
from the collective kinetic energy E]^ij^ for the particle density oscillations. The collective 
kinetic energy is derived as 

1 f 1-2 

-Ekin = -"^Po / dr v"^ = -Bl . (45) 
For the compression modes L = 0, the mass coefficient Bq is given by Q] 

Bo = Bo{q) = {l/2)mpoR',Xo' [l - foi^) - D, jf{x)]^^^^ , for L = 0, (46) 
where xq = qRo- 
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V. RESULTS AND DISCUSSIONS 



We have carried out calculations for the ISGMR in the frameworks of HF based RPA 
and the semiclassical Fermi liquid approach as briefly outlined in the preceding sections. We 
have evaluated the smeared FLA strength function S^^^{E) of Eq. ()42|) and the HF-RPA 
smeared strength function S{E) of Eq. (fT^ for the ISGMR in several nuclei, for /^(r) from 
Eq. (j2ni)- In the subsequent discussions, the HF-RPA results presented are obtained using 
the Skyrme force SkM* I2I] . For the RPA calculations to be highly accurate we discretized 
the continuum in a large box of size 90 fm and use a smearing parameter 7 = r/2 = 1.0 MeV, 
in evaluating the RPA Green's function (see Eq. ()15|) ), and allow particle- hole excitations 
up to 500 MeV (see Ref. jOl for the details). In case of the FLA calculations we have 
adopted the values of po = 0.14 fm'^, ep = 32.85 MeV and Rq = 1.2 ■ A^^^ fm. In Table HI 
we compare the FLA and RPA results for the centroid energies Eqi and E02 corresponding 
to the main and overtone mode of the ISGMR, respectively. We see from this table that the 
FLA and RPA results are in qualitative agreement. The small differences (< 10%) can be 
understood by the fact that the centroid energy mainly depends on the size of the system. 
The values of the mean square radii and the higher moments of the ground state density 
distribution are smaller in the FLA than the ones obtained from HF calculations. Note also 
that the ratio Eq2/Eqi in both models considered is greater than two (~ 2.2 — 2.4). We 
point out that E02 > 2i?oi is due to the Fermi-surface distortion effect as noted earlier in 
Ref. 0. 

We have performed a comparison of the macroscopic transition density p™'^'^'^(r) with the 
ones obtained within the HF-RPA, p^^^^^^(r), and the FLA, p^^'^{r), approaches for the 
main resonance L = and its overtone. The FLA transition density is given by Eq. ()33|) 
with the wave number q obtained from the secular equation (j36|) . The contribution of the 
ISGMR overtone to the EWSR for the case of the transition operator /^(r) is given by 

mo2(0 = E02 J dr fdr)pl%ir) 

where £'02 = huj{qo2)- The eigenfrequency c<j(go2) is obtained from the dispersion equation 
and the wave number go2 is the second (overtone) solution to the secular equation (jSHI)- 
The EWSR for the transition operator /g(r) reads 
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(47) 



Now, we will determine the parameter ^ from the condition that the value of the ratio of 
the partial sum mo2(0 to the total sum mi(^) is the maximum. In Fig. 1, we have plotted 
the ratio ^02(O/''^i(O a function of the parameter ^ for the nucleus ^°'^Pb obtained in 
the FLA (dashed line) and the HF-RPA (solid line) models . As seen from Fig. 1, the 
maximum value of the ratio of ^«02(0 to «^i(^) for the FLA model is achieved for ^ = 68.3 
fm^ where the overtone exhausts about 73% of the EWSR. In the case of HF-RPA, the 
maximum ratio is achieved for ^ = 78.6 fm^ where the overtone, in the energy range of 
25 - 50 MeV, exhausts about 60% of the EWSR. It is interesting to note that if we use 
the FLA and HF ground-state densities to calculate ^ from Eq. (jH^ . we get ^ = 72.2 
and 79.0 fm^, respectively. These values are close to the corresponding ones (68.3 and 78.6 
fm^) obtained by the condition of maximizing the ratio ?t2,02(O/"^i(O- This means that the 
mixing parameter ^ in the transition operator /^(r) of Eq. (f^^ can also be derived from 
the condition that the main ISGMR gives a minimal contribution to the energy-weighted 
sum rule mi(^). The FLA as well as HF-RPA calculations show that the difference in the 
values of ^ obtained in both conditions does not exceed ~ 0.2%. We have also calculated 
the dependence of the parameter ^ on the nuclear mass number A. Following the same 
procedure we find for ^°Zr, ^^^Sn and ^^^Sm nuclei the value of mixing parameter ^ = 38.6 
(48.5), 45.9 (56.7) and 53.1 (66.1) fm^ from the FLA (RPA) calculations, respectively. These 
values can be well approximated by ^ = 1.89A^/'^ and 2.36A^/^ fm^ for the FLA and RPA 
approaches, respectively. 

In Fig. 2, we plot the FLA and RPA results for the fraction energy- weighted transition 
strength as a function of the excitation energy obtained for the transition operator /^(r) 
for the ^°^Pb nucleus. We use ^ = 68.3 and 78.6 fm^ in the FLA and RPA calculations, 
respectively. One can clearly see that the RPA calculation yields a wide resonance of width 
of ~ 10 MeV around the excitation energy 30 — 35 MeV which corresponds to the ISGMR 
overtone. Whereas, in the case of FLA, the transition operator /^(r) (Eq. (j29|) ) with an 
appropriate value of ^ gives rise to a well defined resonance for the overtone mode. We 
also notice that the RPA results have the reminiscence of the ISGMR main tone but it is 
practically eliminated in the FLA calculations. 

In Figs. 3a and 36, we compare the radial macroscopic transition density p'^'^'^^{r) of 
Eq. (pH) . obtained using the HF ground-state density, and the corresponding FLA and HF- 
RPA ones for the overtone of the ISGMR. The radial transition density ptrif) for a certain 
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multipolarity L is given by ptr(i") = Ptr('")^Lo(f). The macroscopic transition densities for 
the overtone in Figs. 3a and 3b are not the same. Because, for an appropriate comparison, 
in Fig. 3a, we have plotted p™'^'^{r) obtained using ^ = 78.6 fm^ and it is normahzed to 60% 
of the EWSR. On the other hand, P^^'^^{t) plotted in Fig. 3b corresponds to ^ = 68.3 fm^ 
and is normalized to 73% of the EWSR. Further, the RPA transition density is calculated 
by averaging over the energy range of 25 - 50 MeV. Notice that the shift of the nodes of 
Ptr^02(^) ^'^ the left with respect to the ones of p^^'^^ij') is caused by the fact that in contrast 
to Eq. ()30|) used for the macroscopic transition density, in the FLA we use sharp nuclear 
surface, see Eq. (jHSI)- We also looked into the energy dependence of the RPA transition 
density for the operator /^(r) over the range of energy employed for the averaging. We see 
that over the entire range considered, the transition density has two nodal structure and the 
distance between the nodes decreases with the increase in energy. For example, averaging 
over AE = 0.5 MeV range, we find that at the excitation energies of 30, 40 and 50 MeV, 
the distances between the two nodes are 3.1, 2.7, and 2.4 fm, respectively, which reflects the 
fact that the microscopic transition density is state dependent. 

We have used the microscopic transition densities for the operator /^(r) to evaluate the 
cross-section for the ISGMR overtone mode excited via inelastic scattering of a-particles 
with energies 240 and 400 MeV. We have used the folding model (FM)-DWBA to calculate 
the excitation cross-section (see Ref. 



23| for details). We find that for the a-particles 
with 400 MeV energy, the calculated cross-section is about 7-10 times higher than the one 
obtained for a-particles with 240 MeV. Note that for a monopole resonance the cross-section 
is maximum at 0°. The values of the cross-section at 0° for the peak energy of the ISGMR 
overtone are 0.5 and 3.5 mb/(sr MeV) for the case of 240 and 400 MeV, respectively. We 
point out that the maximum cross-section for the case of 240 MeV a-particles is below the 

n 

current experimental sensitivity of about 2 mb/(sr MeV) 24]. It may be possible to identify 
the ISGMR overtone mode with 400 MeV a-particles. 

The transition density ptr(r) for the compression modes is distributed over the nuclear 
interior and has a node close to the nuclear surface for both the main ISGMR and its 
overtone. The transition density Ptr,02(r) of the overtone has an additional node in the 
nuclear interior. This feature of Ptr(r) can be tested by evaluating the strength distribution 
So{k) of the electromagnetic operator jL{kr)YLo{f). The strength function Son{k) for a 
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certain eigenstate n is given by 

Sonik) = \Ion{k)f (49) 

where 

hnik) = J dr Ptr,on(r) jo{kr)Yoo{f). (50) 

The strength function S'o„ is related to the excitation function of electron-nucleus scattering 
in the Born approximation. We use Eqs. ()49|) and ()50|) to calculate the energy-weighted sums 
mQi{k) and mQ2{k) for the main ISGMR and its overtone, respectively. In Fig. 4, we display 
the k dependence of the fraction energy- weighted sums moi{k)/mi{k) and mo2{k)/mi{k) for 
the ^ospb nucleus obtained from the FLA (dashed line) and HF-RPA (solid line) approaches. 
It can be seen from Fig. 4 that mQi{k)/mi{k) and mQ2{k)/mi{k) depend strongly on k. A 
shift of the maximum of the ratio mQ2{k) / mi{k) for the overtone to the higher value of wave 
number k is due to the more complicated nodal structure of the transition density associated 
with the overtone as compared with the main resonance. This shift can be exploited to 
separate the modes in electron-nucleus scattering. In Fig. 5 we plot the surface and the 
volume contributions to the integral in Eq. (jKn|l for the transition density associated with 
overtone mode (see Eq. (j33p ). For smaller fc, there is a cancellation between the surface and 
the volume contributions leading to a peak structure for the overtone response as shown in 
Fig. 4. 

VI. SUMMARY AND CONCLUSIONS 

Starting from the local strength function S'(r, E) and using the smearing procedure, 
we have extended the quantum expression for the transition density ptr,n(r) to the case of a 
group of the thin-structure resonances which are localized in the GMR region and are excited 
due to the specifically chosen transition operator F({rj}) = Yld=i fi^i)- approach was 
applied to the study of the transition density of the ISGMR overtone. In this case, an 
appropriate form of the transition operator /(r) is given by /(r) = /^(r) = (r^ — ^r^)yoo(f), 
see Eq. (j29|) . The mixing parameter ^ was determined from the condition that the transition 
operator /^(r) provides for the single overtone the maximum fraction of the energy- weighted 
sum rule mi of Eq. (PT|) . The mixing parameter ^ depends on the nuclear mass number A. 
This dependence is well approximated by ~ 2 ■ A^^^ fm^. 
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We have applied our smearing procedure (using /^(r) associated with the maximum 
FEWSR of the overtone) to the evaluation of the smeared out transition density ptr,R(r) of 
Eq. within the HF-RPA. We have shown that the smearing procedure for the ISGMR 
overtone region provides a simple two nodal structure of ptr,R(r) (see the solid line in Fig. 3a), 
as expected for the L = overtone. Moreover, the transition density pti. R(r), obtained by 
the averaging over many quantum states, resembles its macroscopic counterpart. This fact 
is well illustrated in Fig. 3a by comparing the quantum smeared transition density ptr,R(r) 
with the macroscopic one p^^^^{r) of Eq. (|3(J|). An independent derivation of the smeared 
out transition density ptr,R(r) can be also obtained using the semiclassical approaches. In 
Sec. IV, we have applied a simple semiclassical Fermi-liquid approximation to the evaluation 
of the smeared out (in quantum mechanical sense) transition density pf^^{r). We have used 
the same form of the transition operator /^(r) as in the case of quantum HF-RPA calculation 
to provide an additional check of the derivation of the mixing parameter ^ from the energy- 
weighted sum mi. We found a good agreement between the values of parameter ^ obtained 
in both the quantum and the semiclassical approaches. It is important to emphasize that 
equation ()32|) together with Eq. ()20|) provides a simple expression for the macroscopic 
transition density that can be employed in the folding model-DWBA analysis of excitation 
cross-section of the ISGMR overtone. 

The nodal structure of the semiclassical transition density pf^^{r) is similar to that of 
both the quantum, ptr,R(r), and the macroscopic, p^^'^^{r), cases, see Figs. 3a and 3b. A 
discrepancy occurs in the surface region, where the particular behavior of p^^^{r) is due to 
the assumption of the sharp surface of the nucleus in the FLA model. This discrepancy is not 
so significant in the integral quantities like the strength functions. This is illustrated in Fig. 
4 for the case of the nuclear response to the electromagnetic-like external field ~ jo(^'")^oo(f). 
The ratios mQi{k)/mi{k) and mQ2{k)/mi{k) for the ISGMR and its overtone, respectively, 
show an distinct feature in the fc-dependence. Namely, for a certain value of the wave number 
k, the strength function for the overtone reaches a maximum whereas the contribution of the 
main resonance to the strength function is strongly suppressed. This fact can be exploited 
to separate the ISGMR and the overtone modes in electron-nucleus scattering by varying 
the electron's momentum transfer k. 
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Figure captions 



Fig. 1. The ratio ^7^o2(0/'"^i(0 ^^e partial contribution of the overtone to the EWSR 
as a function of the parameter ^ in the transition operator /^(r) (Eq. ()29p) obtained within 
the FLA (dashed hne) and the HF-RPA (sohd hne) approaches, for the monopole mode 
L = in the nucleus ^°'^Pb. 

Fig. 2. The FLA and HF-RPA results for the fraction energy- weighted transition strength 
for the operator /^(r) with ^ = 78.6 MeVfm^. 

Fig. 3a. The HF-RPA transition density, p^,^~^^^{r), multiplied by Airr'^ for the overtone 
of the ISGMR in the nucleus ^osp]^ (solid line) and the corresponding macroscopic transition 
density p^^^^{r) taken at ^ = 78.6 fm^ (dotted line). 

Fig. 36. The FLA transition density, p^^^{r), multiplied by 47rr^ for the overtone of the 
ISGMR in the nucleus ^°^Pb (dashed line) and the corresponding macroscopic transition 
density p™^'^''(r) taken at ^ = 68.3 fm^ (dotted line). 

Fig. 4. The ratio mQn{k)/mi{k) for the main (n = 1) ISGMR and its overtone (n = 2), 
respectively, as a function of the wave number k obtained for electromagnetic operator 
jo{kr)YQQ{f) for the ^osp]^ nucleus. The dashed and the solid lines represent the FLA and 
HF-RPA results, respectively. 

Fig. 5. Partial contributions of the volume ("vol") and surface ("surf") terms of the FLA 
transition density of the ISGMR overtone (see Eq. (|H3j)) to the integral in Eq. (j^UI) (dashed 
lines). The solid line shows the sum of both the volume and the surface terms. 
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TABLE L Comparison of the FLA and RPA results for the centroid energies (in MeV) for the 
main {Eqi) and the overtone (-E02) modes of ISGMR. In the case of RPA calculations the values 
of Eqi are obtained by integrating the strength function for the operator r^Yoo over the energy 
range of - 60 MeV and the values of £'02 are obtained using the operator (r^ — ^r^)yoo(r) and 
the energy ranges of 35 - 60, 28 - 60, 27 - 55 and 25 - 50 MeV for ^^Zr, ^^^Sn, ^^^Srnand '^°^Pb 
nuclei, respectively. The experimental data for the main tone is taken from the Ref. 



Nucleus 


FLA 


RPA 


EXP. 


Eqi Eq2 Eq2/Eqi 


Eqi Eq2 Eq2/Eqi 


-E'oi 


90 Zr 
ii6Sn 

208pb 


19.6 43.0 2.2 

18.2 40.4 2.2 
17.0 38.5 2.3 

15.3 35.5 2.3 


18.1 43.8 2.4 
16.5 39.1 2.4 

15.7 36.8 2.3 

13.8 33.7 2.4 


17.89±0.20 
16.07±0.12 
15.39±0.28 
14.17±0.28 
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